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Abstract 

Using models ranging from schematic one with a simple quark distribution to more realistic blast 
wave, we study the elliptic flow of hadrons produced from coalescence of quarks and antiquarks 
in the quark-gluon plasma that is formed in ultrarelativistic heavy ion collisions. In particular, 
we study effects due to azimuthal anisotropy in the local transverse momentum distribution of 
quarks, as generated by their position-momentum correlations as a result of radial flow and/or jet 
quenching. We find that even if quarks have large local non-elliptic anisotropic flow, the elliptic 
flow of produced hadrons can still scale with their constituent quark numbers. This scaling is, 
however, broken if the radial flow of coalescing quarks is anisotropic and/or if the momentum 
dispersion of quarks inside hadrons is included. 

PACS numbers: 25.75.-q,25.75.Ld,25.75.Nq 
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I. INTRODUCTION 



In searching for the signatures of the Quark-Gluon Plasma (QGP) produced at the Rel- 
ativists Heavy Ion Collider (RHIC) and to study its properties, the elliptic flow, i.e., the 
second harmonic in the azimuthal distribution of the transverse momenta of produced par- 
ticles, has plaved an important role. It was shown in Ref.jl| using the parton cascade model 
and in Refs. P, [J using the hydrodynamic model that the elliptic flow in heavy ion colli- 
sions at RHIC is sensitive to the properties of formed partonic matter. The observed large 
elliptic flow of charged hadrons at RHIC is indeed consistent with an initial partonic 
matter that is not only dense but also strongly interacting 6]. For identified midrapidity 
hadrons except pions, it was further observed that the dependence of the elliptic flow on 
transverse momentum becomes similar if both are divided by the number of constituent 
quarks in the hadrons, i.e., two for mesons and three for baryons. Such scaling of the el- 
liptic flow of hadrons according to their constituent quark numbers [7 J8 is consistent with 
the picture of constituent or massive quarks coalescing into hadrons 3, 111 U, 12. ll^. In 
particular, the quark number scaling of hadron elliptic flow is exact in the so-called naive 
coalescence model, which only allows quarks and antiquarks (generally called quarks in the 
rest of the paper) with same transverse momentum to coalescence into hadrons. This scaling 
is, however, violated in more realistic quark coalescence models that takes into account the 
momentum distribution of quarks inside hadrons or higher Fock states in 

hadron wave functions jlfij . According to Ref. Q|, the former leads to a violation of the 
quark number scaling of hadron elliptic flow by about 10% for mesons and 15% for baryons, 
resulting in a relative baryon/meson violation of about 5% which is similar to that observed 



in the most recent experimental results |17j]. The hadron elliptic flow is also affected by 
resonance decays 3|. Although this effect has been estimated to be small, including the 
resonance decay effect can largely account for the observed quark number scaling violation 
in the pion elliptic flow HQ. 

Other sources for the violation of the quark number scaling in hadron elliptic flow have 
also been studied recently. For example, it was shown that anisotropics in the quark phase- 
space density and effective emission volume at hadronization can lead to the breaking of 
hadron elliptic flow according to their constituent quark numbers Furthermore, it was 
pointed out in Ref. |2Qj that in usual derivation of the quark number scaling of the elliptic 
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flow based on the naive quark coalescence model 

BUI 111 ] . correlations between positions 
and momenta of quarks, as expected from the hydrodynamic or the parton cascade model, 
were ignored as only the spatially averaged quark momentum distributions were used. Since 
the position-momentum correlations can generate large non-elliptic anisotropy in the quark 
local momentum distribution, which may not be averaged out after integrating over the 
quark spatial distribution, the quark number scaling of hadron elliptic flows may thus be 
strongly violated. 

To understand more clearly what role is played by quark position-momentum correlations 
and how the quark number scaling of hadron elliptic flow appears in the quark coalescence 
model, we examine in the present paper the relation between the elliptic flows of quarks and 
mesons in models that range from a simple schematic model to the more realistic blast wave 
model. Our analyses indicate that although large azimuthal anisotropy generally exists in 
the quark local transverse momentum distributions in these diverse models as well as in the 
model used in Ref.[12j, it does not necessarily lead to the violation of the quark number 
scaling of hadron elliptic flows if the quark global transverse momentum distribution does 
not have odd-order anisotropic flows. 

II. ANISOTROPIC FLOW 

Following the notations of Ref.Q], the spatial and transverse momentum distributions 
S(x,p T ,y) of particles produced in relativistic heavy ion collisions can be expanded in 
Fourier series as 

S(x, pr, y) = S (x, PT , y) jl + 2 £ Re (c n (x,p T , y)e~ m ^) j , (1) 

where px = Pt(cos0 p , sin</> p ) is the transverse momentum of a particle with cf) p denoting its 
azimuthal angle with respect to the reaction plane, and y is its rapidity. The Fourier coef- 
ficient c n (x,pT,y) denotes the nth-order local anisotropic flow of particles with momentum 
Pt in an infinitesimal space-time volume around x. In general, c n is complex and is written 
as c n = v n + iu n , with the real and imaginary parts given by averaging over the azimuthal 
angle 4> p of particle momentum, i.e., v n = (cos(n0 p )) and u n = (sxn(n<p p )) . 

Experimentally, only the spatially averaged global anisotropic flows are measured, i.e., 
r (7>r v) = (r (r ^ v)) - I d ' xc n(^PT,y)S(x,p T ,y) 



The real part of the two lowest orders v\ and v 2 are the directed and elliptic flows, respec- 
tively. The imaginary part u n always vanish in collisions between identical spherical nuclei 
due to reflection symmetry with respect to the reaction plane, i.e., invariant under the 
transformation — > —0. Furthermore, for particles at midrapidity in collisions with equal 
mass nuclei, anisotropic flows of odd orders vanish as a result of the additional symmetry 
— > cf) + 7r. In the present study, we consider only particles at midrapidity in collisions of 
identical spherical nuclei. 

In the naive coalescence model , which only considers the momentum distri- 

bution of quarks, thus neglecting correlations between quark positions and momenta, and 
allows quarks with same transverse momentum to coalesce into hadrons, one obtains a sim- 
ple relation between meson and quark elliptic flow if all v n of quarks vanish except for n = 2; 
i.e, 



Neglecting v\ in the denominators due to its small value, a scaling of the elliptic flow of 
mesons according to their constituent quark number then follows. 

Including quark position-momentum correlations, the meson local elliptic flow v 2) m is 
then related to quark local v n and u n via 



_ 2t> 2 + v\ - u\ + 2 Y^ =l {v k v 2+k + u k u 2+k ) , , . 

V%M 1 + + ' ( ' 



Although higher-order quark local anisotropic flows v n and Un with n > 3 are likely small, 
quark local Vi and u% can be large as first pointed out in Ref.|2(|. It is not obvious that in the 
spatial average the contribution of quark local directed flow t>i, U\ vanishes. The observed 
approximate scaling of hadron elliptic flow according to the constituent quark number may 
simply due to accidental cancellations of quark local directed flow as suggested in Ref.[20j. 
In the present paper, we study the effect of nonvanishing quark local directed flow on the 
meson elliptic flow in various models. Our analysis shows that after spatial average its effect 
either disappears or is small, thus having negligible or only small effect on the elliptic flow 
of hadrons. 
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III. A SCHEMATIC MODEL 



We first consider a schematic model in which quarks are distributed uniformly in space 
with their momenta pointing radially outward, i.e., same azimuthal angles 4> p and 4> r for the 
momentum and position vectors of a quark. This quark distribution is qualitatively similar 
to that of quenched jets produced in heavy ion collisions at RHIC j^l, 22|. Including also an 
azimuthal anisotropy of the form cos(20 p ) in the quark transverse momentum distribution, 
the quark distribution becomes 

dN 

oc (1 + 2a 2 cos(20 p )) 5(4> r - <j) p ), (5) 



d 2 rd 2 p 

where the coefficient a 2 is the quark elliptic flow v 2yq . 

Considering only quarks with spatial azimuthal angles less than <fi' r , the Fourier coefficient 
c n>q of resulting momentum distribution can be evaluated by integrating r from to (f>' r , 
i.e., 

1 r<t>'r r 2n 



N Jo Jo 



d<j) p [1 + 2a 2 cos(20 p )] 5{<f) r - p )e m ^ 

t -1 \ 
r I 

(6) 



ill 2n e m,p r 

T7 < ~ + a 2 — 2a 

A n n z — 4 n 



e i(n-2)4>' r e i(n+2)<P' r 

+ 



n-2 n + 2 

where N = (p' r + a 2 sin(20^.) is the normalization factor shown in the denominator of Eq. 
For n = 2, the above equation gives 



1 

N 



C2, q = — i a 2 (j) r + i 



2 + 4 v2 + 4 ] 



(7) 



The global quark elliptic flow is obtained from Eq.fJZJ) by letting <p' r = 2ti, and this gives 
V2,q = 0,2 as expected. It can be checked that all higher-order global quark anisotropic 
flows vanish as the first two terms in Eq.fjHJ) are exactly canceled by the last two terms 
when (j)' r = 2tt. On the other hand, for quarks in a localized space the distribution contains 
an infinite number of Fourier coefficients. For example, for quarks in the first quadrant, 
which corresponds to integrating (f) r from to (j)' r — 7r/2, all c ra 's are nonzero. Assuming an 
anisotropy parameter a 2 = 0.1, we find that for quarks in the first quadrant, the local directed 
flow is v\ = 0.68 and its imaginary part is u\ = 0.59. Both are close to those from the MPC 
code for Au+Au collisions at yfs = 200 GeV with b = 8 fm j^. Except for the local elliptic 
flow v 2>q = v 2>q = 0.1, other real even Fourier coefficients all vanish, and this is similar to 
that in the midrapidity region of heavy ion collisions at RHIC. The imaginary even Fourier 



coefficients are in general nonzero, and we have for the second Fourier coefficient w 2j(? = 0.64. 
We note that the large azimuthal anisotropy in the quark local momentum distribution in 
the present model is caused by the strong position and momentum correlations in the quark 



distribution, which is also present in the blast wave model 



23| as well as in the transport 



model j^J. These models further include, however, effects due to thermal motions of quarks. 
In the naive coalescence model, the transverse momentum spectrum of produced mesons 



is given by 

dN M 
d 2 p 



12 



9m 



(8) 



' d 2 T i d 2 'p 

where g m is the statistical factor for a quark and an antiquark to form a colorless meson 
with certain spin and isospin. 

With the quark distribution of Eq.fjSJ), carrying out the integral in Eq.fjHJ) over the az- 
imuthal angle from to (f)' r gives following meson anisotropic flow in a restricted space: 
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(9) 



where N 
space is then 



l + 2a 2 ) 



^ + 2a 2 sin(2<^) + a!/2sin(4</£). The local elliptic flow in this restricted 



C2,M 
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Including all quarks, i.e., letting <p' r = 2ir in Eq.Q, cancellations similar to those in 
the evaluation of the quark global anisotropy flows lead to the following meson transverse 
momentum distribution: 

dN M 



d 2 p 



oc 1 + 2a 2 + 4a 2 cos(20 p ) + 2a 2 cos 



The above equation shows that the meson elliptic flow is 

2lJ 2,g 



V2,M 



1 + 2^' 



and there is also a 4th-order anisotropic flow given by 

V4,M 



2,9 



1 + 24/ 



(11) 



(12) 



(13) 



These results are exactly what one expects from the naive coalescence model when all quark 
anisotropic flows c n except v 2 are zero [n| . 

The above result holds for any value of the anisotropy parameter a 2 = v 2)Q , and thus for 
any value of local anisotropic flow c n>q . In other words, despite the fact that locally there 
can be very large anisotropy in the quark momentum distribution, they are not relevant for 
the global elliptic flow of mesons that are produced from the coalescence process. The quark 
number scaling of hadron elliptic flow is exact as expected when the quark distribution has 
only non-zero global elliptic flow. 



IV. THE BLAST WAVE MODEL 



The schematic model introduced above can be made more realistic by including thermal 
motions of quarks. Adding a random component to the quark momentum weakens the 
extreme angular correlation between the positions and momenta of quarks in the schematic 
model. As shown in the following, mesons produced from the naive coalescence remain to 
have an elliptic flow that follows the quark number scaling if the radial flow is isotropic, 
i.e. independent of azimuthal angle. The large anisotropy in the quark local momentum 
distribution, particularly the V\ tq and Ui >q , does not affect the scaling relation between quark 
and meson elliptic flows. On the other hand, the quark number scaling of hadron elliptic 
flow can be violated when the quark radial flow is anisotropic, i.e., its velocity varies with 
the spatial azimuthal angle <p r . 



A. Blast wave with isotropic radial flow 



We first consider the case that quarks have a relativistic Boltzmann distribution with 

n 

isotropic radial flow velocity [25||, i.e., 



dN, 



oc exp 



~r coshp + ^ sinhpo cos( 



[1 + 2a 2 (p T ) cos(20 p )] , (14) 



d 2 p T d 2 r 

where = (p^+m 2 ) 1 ^ 2 . In the above, dependence on the spatial and momentum rapidities 
y and rj has been dropped as we consider a Bjorken boost invariant distribution at y — r] — 0. 
For simplicity, both the quark density distribution and their flow profile (3t = tanh po 
are taken to be independent of the radial position. Also, the spatial shape of the quark 
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distribution is assumed to be cylindrical. Although spatial ellipticity can cause the breaking 
of quark number scaling of hadron elliptic flow as discussed in Ref.[l9|, it is not essential for 
the purpose of present study. 

For quarks with spatial azimuthal angle less than <j)' r , the different harmonics in the 
Fourier expansion of the quark distribution in Eq. (jl4j) are given by 

1 ^2tt 

N Jo 'JO 

x[l + 2a 2 (p r )cos(20 p )], (15) 



1 f27T 

c n ,q{PT) = -t-j d<p p d(/) r exp [a(p T ) cos(0 r - <p p ) + inc, 
JS Jo Jo 



with N = 2ir{I [a(pT)]4>' r + a2(pT)h[ot(pT)] sin(20j,)} and a(pr) = (pt/T) sinhp - If ah 
quarks are included, i.e., <p' r = 2tt, then all orders of anisotropic flow except V2 vanish. This 
can be easily seen as the spatial integral can now be expressed in terms of the Bessel function 
Iq and is thus canceled by the same factor in the normalization factor N in the denominator. 
What is left is then the integral over <p p which is non-zero only for n = 2. 
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FIG. 1: Quark local anisotropic flow as a function of transverse momentum from averaging over 
quarks in the first quadrant of the transverse plane. 



On the other hand, because of quark radial flow, none of the harmonics in the Fourier 
expansion vanishes if only quarks in a limited space are considered. To estimate the values 
of the different harmonics in the present model for heavy ion collisions at RHIC, we take 
a slope parameter T = 170 MeV for the quark transverse momentum distribution, which 
is consistent with the phase transition temperature from lattice QCD calculations. Masses 



of quarks are taken to be m q 



m- n 



300 MeV as in Ref. 



12 1 . For the strength of quark 



anisotropic flow, we take po — 0.4 and 02 (pt) = 0.1 tanh(1.2pr) i n order to have radial and 
elliptic flows with magnitude similar to those observed in experiments. Results for the first 
two local harmonics Vi q , ui jq , v 2a and u 2%q in the first quadrant are shown in Fig IJ and they 
are seen to have values similar to those from the MPC quark cascade model 24 1. 

The meson spectrum obtained from the naive coalescence of quarks, i.e., only quarks in 
the same phase space point can coalesce, can be easily derived from the quark distribution 
of Eq. (fT4*|) by integrating over both spatial azimuthal angle (f) r and momentum azimuthal 
angle p , and the result is given by 



oc e - 2&(pT/2) I [2a(pr/2)] 



dN M 
d 2 p T 

x 1 + 2a 2 2 (p T ) + 4a 2 (pr) cos(20 p ) + 2a 2 2 {p T ) cos(40 p )] , (16) 

where = (m T /T) cosh p . Comparing to Eq.(JTTJ), we see that the quark number scaling 

of meson elliptic flow appears as in the schematic model of Section 11111 This example thus 
demonstrates again that the naive coalescence model leads to the quark number scaling 
of hadron elliptic flow, independent of the presence of large anisotropy in the quark local 
momentum distribution. 

If we had considered quarks in a finite rapidity range, integrating over the rapidity leads 
to the Bessel function Ki(a) instead of the e~ a in Eq.([TBJ). This does not, however, affect 
the conclusion from the simpler picture considered in the above. 



B. Blast wave with anisotropic radial flow 



The above blast wave model can be made more general by allowing the radial flow 
velocity to depend on both the azimuthal angle in space and the transverse momen- 



tum. Specifically, we introduce a radial flow profile /3;r(pT,0r) = tanh p(pr, 



with 



p(pT, 4> r ) — Po + P 2 (pt) cos(20 r ). This is similar to that in the blast wave model of Ref.| 
except for the pt dependence of the azimuthal anisotropy parameter p 2 . The latter makes 
it possible to parameterize the typical increase and then saturation (or slight decrease) of 
the elliptic flow with transverse momentum, that is observed in experiments jj, ^| and pre- 



dicted by transport models 



27 



28 



29j | . With anisotropic radial flow, quarks can acquire an 



elliptic flow without introducing explicitly anisotropic azimuthal distributions in space and 
momentum as in the models discussed so far. The quark distribution in the present model 
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is then given by 

(17) 



dN qA 

oc exp 



d 2 p T d 2 r 



'' lT COsh p(p T} <f> r ) + Sinn P(PT> 0r) COs( 



Defining a{jpT , 4>r) = ( m r/T) cosh p(p T , <j> r ) and a(pT,4>r) — (Pt/T) sinhp(p T , r ), the 
transverse momentum dependence of the quark elliptic flow can be easily calculated as in 



Ref. 26], and the result is given by 

- / \ Io W #rexp[-«(p T ,0 r .)]J 2 [a(p T ,0 r )]cos(20 r ) 
J d<p r exp [-a(p T , r )J I [a{p T , (p r )\ 

In the naive quark coalescence model, which only allows quarks with same momentum 
to coalesce into hadrons, the meson distribution is simply given by the following one- 
dimensional integral: 

dN M r 2n 

—p— oc d(j) r exp [-2a(p T /2, <p r )] exp [2a(p T /2, <p r ) cos(<p r - <p p )] . (19) 

The elliptic flow of mesons is then 

d<j> r exp [-2a(p T /2, r )] I 2 [2a(p T /2, <f> r )] cos(20 r ) 
2 ' MlPTj / o 27r C /0 r exp[-2«(p T /2,0 r )]J o [2a(p r /2,0 r )] ' 1 J 

Contrary to the case of isotropic radial flow, the quark number scaling of hadron elliptic 

flow is not obvious when the radial flow is anisotropic, and a violation is generally found as 

shown below. 

Using the anisotropic function p 2 = 0.25/{l + expQ>r — 0.5)/0.75]} + 0.014 in the quark 
radial flow, we have evaluated numerically both the quark and scaled meson elliptic flows, 
and they are shown in Fig|21 by the solid and dashed thick lines, respectively. It is seen 
that the quark number scaling of hadron elliptic flow is appreciably broken (about 20%) at 
Px/n ~ 1.5 GeV. Furthermore, the meson elliptic flow in the present model is larger than 
the one expected from the quark number scaling, in contrast with the smaller meson elliptic 
flow from most other sources of quark number scaling violation 15, 16, 19]. Since the 



local harmonics in the quark momentum distribution in the present case have values very 
similar to those in the previous blast wave model with isotropic radial flow, they are not 
the reason that the quark number scaling of meson flow is violated. The violation is caused 
instead by the larger coalescence probability for quarks in the reaction plane than out of the 
reaction plane as a result of anisotropy in the radial flow. This can be seen by examining 
the second harmonic in the spatial distribution s 2 = (cos(20 r )) of quarks and mesons. As 
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FIG. 2: Scaled momentum elliptic flow U2 and spatial elliptic deformation S2 as a function of scaled 
transverse momentum. Solid lines are the quark V2 >q (thick) and «2,g (thin), while dashed lines are 
V2,m (thick) and ~S2,M (thin) of mesons. 

shown in Fig|21 the scaled spatial elliptic deformation S2 of mesons (thin dashed line) is 
indeed larger than that of the quarks (thin solid line), except at low transverse momenta. 

Our results at low transverse momenta needs, however, to be interpreted with caution 
as the coalescence model violates unitary if applied to the entire phase space. Although 
the coalescence model is applicable for describing the production of rare particles such 
as those with intermediate transverse momenta, a more dynamical approach that takes 
into account properly energy and entropy conservation is needed for more abundant low 
transverse momenta particles. 

We would like to point out that for quarks with transverse momentum pr/n > 1.5 GeV, 
their elliptic flow is not likely affected by anisotropy in their radial flow. The momentum 
distribution of these quarks is thus more similar to that in the schematic model of Sect. IHII 
for which the quark number scaling of hadron elliptic flow holds exactly. 



V. HADRON WAVE FUNCTION EFFECT ON ELLIPTIC FLOW 

Besides anisotropic radial flow, the quark number scaling of hadron elliptic flow can also 
be violated by the coalescence of quarks with different momenta and at different positions 
when the hadron wave function is taken into consideration as in Refs.QQ. To see this 
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explicitly, we generalize Eq.© to three dimensional space and replace the resulting two 



delta functions (27r) 3 5^(ri — ^(^(pi — P2) by the meson Wigner function from Ref. 
i.e., 
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(21) 



where A x = A 1 is the width of the quark momentum distribution in the meson. 
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FIG. 3: Scaled elliptic flow as a function of scaled transverse momentum. The solid line is the 
quark elliptic flow V2g- Meson elliptic flow from the coalescence model is shown by open circles for 
a width of hadron Wigner function A p = 0.48 GeV, by filled squares for A p = 0.24 GeV, and by 
filled circles for A p = 0.1 GeV. Dashed lines are drawn to guide the eye. 



Using the same quark distribution function of Ref. [12], which is essentially the blast wave 
model with isotropic radial flow studied in Sect II V M we show in Fig. Elthe result for mesons 
from the coalescence model based on a Monte Carlo evaluation of the coalescence integral. 
Three different values of A,, are used. The filled squares are the result with A p = 0.24 GeV 
that was used in Ref. [12I llq to fit the experimental meson spectra at RHIC. In this case, the 
quark number scaling is violated by about 10% in the scaled momentum region pr/n — 1 — 2 
GeV, as already found in a previous work [12I I15I ] and mentioned in the Introduction. In 



Ref. 



24j . a larger violation of about 20% has been found in the quark number scaling of 



hadron elliptic flow. This may be largely due to the fact that the corresponding value for 
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Ap used in Ref. 24] is about a factor of two larger than our value due to the absence of the 
factor 97r/2 in the hadron Wigner function. Indeed, using A p = 0.48 GeV in the present 
calculation leads to a much smaller scaled meson elliptic flow as shown by open circles, thus 
a much stronger violation of the quark number scaling of the meson elliptic flow. In the 
scaled momentum region pr/n ~ 1 — 2 GeV, the violation of the quark number scaling is 
consistent with the 20% effect seen in Ref.|24|. Also shown in Fig. El by filled circles are 
results obtained with a smaller A p = 0.1 MeV. It is seen that they are now very close to the 
quark elliptic flow given by the solid line, as expected in the limit of A p — > 0. 
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FIG. 4: Local harmonics of quark momentum distributions as a function of transverse momentum 
from averaging over the first quadrant of the transverse plane. 



We have also evaluated the local harmonics of the momentum distribution of quarks in 
the first quadrant. They are shown in Fig. H]for vi )Q (filled squares), Ui >q (open squares), 
V2, q (filled circles), and U2. q (open circles). Both the real and imaginary local directed flows 
Vi >q and U\ A are large, as in the anisotropic blast wave model of Sect II V B*1 shown in FigO 
Our results thus demonstrate that violation of the quark number scaling of hadron elliptic 



flow in the realistic model of Ref. 



12 



3 



is not due to the large local anisotropic flow but is 



rather a result of the finite dispersion of quark momentum inside hadrons. 
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VI. SUMMARY 



To understand the observed quark number scaling in the elliptic flow of identified hadrons 
at RHIC, we have studied a number of models that range from the schematic one with 
a simple quark distribution to the realistic blast wave model that is consistent with the 
underlying hydrodynamic and transport models. In the naive coalescence model, in which 
only quarks with same momentum can coalescence, all models, except the blast wave model 
with an anisotropic radial flow, lead to scaling of hadron elliptic flow according to the 
constituent quark number. In the anisotropic blast wave model, the violation of this scaling 
is due to the azimuthal asymmetry in the quark coalescence probability as a result of the 
azimuthal angle dependence of the radial flow. In all cases, large local anisotropic flows are 
present, but they do not play any role in the final hadron elliptic flow. We have further 
demonstrated in a blast wave model that the quark number scaling can also be violated in 
a more realistic coalescence model which takes into account the momentum dispersion of 
quarks in hadron wave functions. Since effects due to anisotropic radial flow and hadron 
wave functions are small if realistic values are used for these quantities, observation of the 
quark number scaling of hadron elliptic flow in experiments thus provides a unique signature 
for hadronization via the quark coalescence as well as for the existence of the quark-gluon 
plasma prior the production of hadrons. 
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